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Multi-subband effect on spin precession and spin dephasing in n-type GaAs quantum wells is 
studied with electron-electron and electron-phonon scattering explicitly included. The effects of 
temperature, well width and applied electric field (in hot-electron regime) on the spin kinetics 
are thoroughly investigated. It is shown that due to the strong inter-subband scattering, the spin 
procession and the spin dephasing rate of electrons in different subbands are almost identical despite 
the large difference in the D'yakonov-Perel' (DP) terms of different subbands. It is also shown that 
for quantum wells with small well width at temperatures where only the lowest subband is occupied, 
the spin dephasing time increases with the temperature as well as the applied in-plane electric field 
until the contribution from the second subband is no longer negligible. For wide quantum wells the 
spin dephasing time tends to decrease with the temperature and the electric field. 

PACS numbers: 72.25.Rb, 72.20.Ht, 71.10.-w, 67.57.Lm, 73.61.Ey 
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Manipulating and transporting electron spins belong 
to the subjects of intense research in semiconductor 
spintronics,"'^'^ which aims to incorporate the spin degree 
of freedom into the traditional electronic devices. To date 
most of the technical proposals to manipulate the elec- 
tron spin are either through an applied magnetic field or 
via a gate voltage that changes the spin-orbital coupling 
in confined semiconductor structures through the Rashba 
effect. "^'^ Recently it has been realized that a strong in- 
plane electric field also provides another path to manipu- 
late the spin dephasing in the presence of the spin-orbital 
coupling. Rashba and Efros showed that an in-plane ac 
electric field can manipulate electron spins efficiently.^ 
Pramanik et al. investigated the effect of an in-plane 
electric field on the spin dephasing in GaAs quantum 
wire through Monte-Carlo simulation and revealed that 
the presence of the strong electric field enhances the spin 
dephasing. ^'^ A comprehensive many-body investigation 
of the hot-electron effect on spin precession and spin de- 
phasing due to the strong in-plane electric field in n- 
type GaAs quantum wells (QW's) has been reported 
recently* in the electric-field and temperature regime 
where only the lowest subband is occupied. In that inves- 
tigation all scattering, ie., the electron-electron, electron- 
phonon and electron-impurity scattering is explicitly in- 
cluded and is calculated self-consistently. It is discov- 
ered however that the spin dephasing rate decreases with 
the applied electric field. In the narrow QW, it is also 
discovered that the spin dephasing rate decreases with 
temperature. ^'^° 

The spin precession and spin dephasing in wide QW's 
should show quite different temperature and electric-field 
dependence as electrons populate more than one subband 
and therefore experience quite different spin-orbital cou- 
pling strength. In the present paper, we study the multi- 
subband effect on the spin precession/dephasing in n- 
type GaAs QW's, where the spin dephasing mainly comes 



from the D'yakonov-Perel' (DP) mechanism, ^^ with and 
without the high in-plane electric field. The inclusion 
of multi subbands allows us to investigate QW's with 
wide well width and also allows us to study the regime 
of higher electric field where the single subband model 
cannot deal with due to the "runaway effect". *'^^ We re- 
veal the intra- and inter-subband scattering to the spin 
dephasing times of each subband. The model of our in- 
vestigation is composed of (100) GaAs QW's of width 
a with its growth direction along the z-axis. An uni- 
form electric field E and a moderate magnetic field B 
are applied along the x-axis (Vogit configuration). Due 
to the confinement of the QW, the momentum states 
along the z-axis are quantized. Therefore the electron 
states are characterized by a subband index n and a two- 
dimensional wavevector k = (kx,ky), together with a 
spin index a. With the DP term (DPT) included, the 
Hamiltonian of the electron in the QW reads: 

H = ^ {(e„k - eE • R)dnn'Saa' + [ff/isB + h„„'(k)] 

-^Ulk^Cn'Ua' +Hl. (1) 



(Tr, 



Here e„k — k-^/2m* -|- (fc^)„/m* is the energy spectrum 
of the electron with momentum k and effective mass m* 
in the n-th subband. cr are the Pauli matrices. R = 
{x,y) represents the position. h„(k) is the DPT which 
serves as an effective magnetic field with its magnitude 
and direction depending on k. It is composed of the 
Dresselhaus term^^ and the Rashba term."^'^ For GaAs 
QW, the leading term is the Dresselhaus one which can 
be written as: 



,x(k) 



jkx{ky - {n\k'^\n))d„ 
■jky{{n\kl\n) - kl)S„ 



,(k) = ^{n\k.\n')ikl-kl) 
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In the above, (n|fc^|n) represents the average of the oper- 
ator — (^)^ over the electronic state of the n-th subband 
and is therefore n^(7r/a)^ under the infinite- well-depth 
assumption and {n\kz\n') is the matrix element of —id/dz 
between n- and n'-th subbands. 7 here is the spin- 
orbital coupling constant.^* The interaction Hamiltonian 
Hi is composed of the Coulomb interaction H^e, the 
electron-phonon scattering Hph, as well as the electron- 
impurity scattering Hi. Their expressions can be found 
in textbooks. ^^'^^ 

In order to study the hot-electron effect on spin de- 
phasing, we limit our system to a spacial homogeneous 
one in order to avoid the additional complicity such as 
charge/spin diffusion. The kinetic Bloch equations in 
such a system are constructed using the nonequilibrium 
Green function method with the gradient expansion^® 
and can be written as: 

PnU-.aa' *=--'-'' V k/^nk,(j(j' /^nk, (jo"' |colii/^nk, (T(t' Iscatt 7 \'^) 

where Pn]s.aa' represent the single particle density ma- 
trix elements. The diagonal terms describe the electron 
distribution functions Pnk.o-o- = fnka- The off-diagonal 
elements /9„]^ j._j. =p*t_ii = Pnk stand for the inter- 

spin-band polarizations (spin coherence). ^^ The second 
terms in the kinetic equations describe the momentum 
and energy input from the electric field E. PnWaa' \coh 
on the right hand side of the equations describe the co- 
herent spin precession around the applied magnetic field 
B, the effective magnetic field h(k) from the DPT as 
well as the effective magnetic field from the electron- 
electron Coulomb interaction in the Hartree-Fock order. 
Pnkcrcr'Iscatt dcuote the electron-impurity, the electron- 
phonons, as well as the electron-electron Coulomb scat- 
tering. The electron-electron Coulomb scattering comes 
from the Coulomb interaction beyond the lowest order, 
i.e., the Hartree-Fock order. It is noted that in writ- 
ing the kinetic equations, we have neglected the inter- 
subband coherence Pnn'kaa'in ^ n') as these terms are 
much smaller than the intra-subband spin coherence due 
to the small "pumping" terms hnn',z{^)- 

It is seen that all the unknowns to be solved appear in 
the coherent and the scattering terms nonlinearly. There- 
fore the kinetic Bloch equations have to be solved self- 
consistently to obtain the electron distribution and the 
spin coherence. By numerically solving the kinetic Bloch 
equations in the self-consistent fashion, one is able to 
obtain the temporal evolutions of the electron distribu- 
tion functions fnv,a{t) and the spin coherence pnu.{t). 
Once these quantities are obtained, all the quantities 
such as electron mobility /i„, hot-electron temperature 
Te.n as well as the spin dephasing rate for electrons in 
each subband can be deduced. For subband n, the mo- 
bility is given by /x„ = X^k^ /"k^ kx/[m*Ej2i^afni^'^]- 
the electron temperature is obtained by fitting the Boltz- 
mann tail of the electron distribution function; whereas 
the spin dephasing rate is determined by the slope of 
the envelope of the incoherently summed spin coherence 
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FIG. 1: Temporal evolution of the spin signal with width a — 
17.8 nm and electron density Ne^4x lO" cm~^ for T = 200 
K and B — 4 T under zero applied electric field. The initial 
spin polarization P = 2.5 %. Upper panel: The incoherently- 
summed spin coherence p„ of n-th subband; Lower panel: 
The electron densities of n-th subband with spin a. The solid 
curves are for the first subband and the dashed ones are for the 
second subband. Note that the scales of the second subband 
are on the right-hand side of the figure. 



The initial conditions at i = are taken to be pnkiO) = 
and electron distribution functions are chosen to be 
those in the steady state under the electric field but 
without the magnetic field and the DPT.^ Specifically 
/nk,<T(0) is the solution of the kinetic equations (3) in 
the steady state, with the spin coherence p„k7 the mag- 
netic field and the DPT set to zero. The implementation 
schemes of the numerical solution of the kinetic equations 
and the initial condition can be found in Ref. 8. The to- 
tal electron density Ne , the applied magnetic field B and 
the initial spin polarization P = X]nk(-/'nki ^ /nk-i)/-^e 
are taken to be 4 x 10^^ cm^^, 4 T and 2.5 % respectively 
throughout the paper. 

We first study the spin dephasing of electrons in a 
GaAs QW with a ^ 17.8 nm. In this QW, the elec- 
trons mainly distribute in the first two subbands under 
the applied electric field up to 3 kV/cm. Therefore, we 
only consider electrons in the first two subbands with n 
being 1 and 2. 

We first focus on the temporal evolution of the spin sig- 
nals. In Fig. 1 we plot the densities of electrons with spin- 
up and -down as well as the corresponding incoherently- 
summed spin coherence of each subband versus time t 
without applied electric field. The background tempera- 
ture T = 200 K. It is seen from the figure that the elec- 
trons in both subbands undergo damped oscillations with 
a same precession frequency and damping rate. Never- 
theless, as (2|fc2|2) = 4(1|A:2|I), from Eq. (2) electrons 



with a given wavevector in the second subband experi- 
ence an effective magnetic field from the DPT with its 
magnitude three times larger than that of the first sub- 
band. At first glance one may expect that the spin dc- 
phasing rate of electrons in the second subband should 
be about 9 times faster than those of electrons in the first 
subband. It is of particular interest to see that the decay 
rates of the spin signals of the two subbands are almost 
identical. 



0.012 



0.008 



:3 



0.004 




100 



t (ps) 



FIG. 2: Temporal evolution of the incoherently-summed spin 
coherence of each subband pn for the same situation in Fig. 1. 
The blue curves are Pn{t) (n = 1 and 2) from the full kinetic 
equations (Noted that pi{t) and P2{t) match each other); The 
green curves are Pnit) without the inter-subband Coulomb 
scattering; The red curves are Pn{t) without the inter-subband 
scattering. Solid curves: pi(i); Dashed curves: p2{t). 

In order to reveal the physics of the identicalness of the 
spin precession and the dephasing of electrons in different 
subbands, we further plot the temporal evolution of the 
spin coherence of the above case in Fig. 2. For compari- 
son we also plot p„ (t) without the inter-subband electron- 
electron Coulomb scattering (green curves) and those 
without the inter-subband electron-electron and electron- 
phonon scattering, i.e., without any inter-subband scat- 
tering (red curves). It is seen from the figure that if 
we count for all the scattering, the spin coherence of 
the different subbands decays identically. However, once 
the inter-subband of the Coulomb scattering is removed, 
the spin coherence of the lowest subband decays a lit- 
tle slower and that of the second subband decays a lit- 
tle faster than that with all the scattering included. If 
one further removes the inter-subband electron-phonon 
scattering, the spin coherence of the second subband de- 
cays very fast and disappears in the first 20 ps while 
pi lasts much longer than the case with all the inter- 
subband scattering included. It is clear to see that as the 
inter-subband scattering is reduced, the identicalness of 
the spin dephasing of electrons in different subband is 



removed. Once the inter-subband scattering is totally 
removed, the spin coherence differs greatly between the 
different subbands. Therefore, the identicalness is due to 
the strong inter-subband scattering. 

To further elucidate the effect of the inter-subband 
scattering on the spin dephasing, we adopt a much sim- 
plified model: The system is simply described by the elec- 
tron number iVi, N2 and the magnetic momentum Mi, 
M2 of the two subbands. In the non-degenerate case, 
the kinetics of the system are: dNi/dt = X]7=i 2 ^ij^j 
and d{N,M,)/dt = Y. 3=1.2 ^^^3(^3^3) " ^^Af./r/ {i = 1, 
2), with aij standing for the inter-subband scattering 
rate and 1/t/ for the spin dephasing rate of the z-th 
subband. It is easily seen that in the absence of the 
inter-subband scattering, the spin signals in the different 
subbands decay with their own dephasing rate. How- 
ever, once a strong inter-subband scattering is presented 
and the electron populations approach to the equilibrium, 
one is able to obtain that an = —Xj/t and a^ = Xi/r 
{j 7^ i) from the equilibrium condition dNi/dt = 
and the detailed balance condition aijNj = ajiNi with 
Xi = Ni/{Ni + N2). 1/t is the inter-subband relaxation 
rate. Substituting these relations into the kinetic equa- 
tions for the magnetic moments Mi, one gets: 



dMi/dt = -X2{Mi-M2)/t-Mi/tI, 

divh/dt = xi{Mi-M2)/t~M2ItI 



(4) 



In the strong inter-subband scattering limit, i.e., 1/t » 
1/t/, the difference of magnetic momentum Mi — M2 de- 
cays with the rate 1/t. That is the magnetic momentums 
Mi of different subbands become almost identical to each 
other in the time scale of t. Taking this fact into account, 
one easily gets the equation that controls the magnetic 
momentum M = Mi = M2: dM/dt = -J2ii^J'^!)M- 
Thus in the strong inter-subband scattering limit, the 
magnetic momentums of different subbands decay with 
the same rate I/t^ = J2i^i/'''i- Therefore, it is under- 
stood that, in the presence of the strong inter-subband 
scattering, electrons hop frequently among the subbands 
and fast exchange the spin signal. Consequently elec- 
trons experience an average DP effective magnetic field 
of the different subbands during the spin precession and 
acquire the same spin dephasing rate. Moreover, we also 
discovered that in the case of strong electric field where 
the DPT can also change the spin precession rates,® a 
strong inter-subband scattering also causes the identi- 
calness of the spin precession rates of different subband, 
regardless from what expected from the DPT analysis 
that the change of the spin precession rate of the second 
subband should be 3 times larger than that of the first 
one. 

We then turn to the temperature dependence of the 
spin dephasing of electrons in QW's with different well 
width. The spin dephasing time (SDT) versus the back- 
ground temperature without applied electric field is pre- 
sented in Fig. 3 for QW's with a = 17.8 nm and a = 12.6 
nm respectively. In the figure we also plot the SDT in 
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FIG. 3: SDT Ts vs. the background temperature T without 
the apphed electric field for two QW's with width a — 17.8 nm 
(•) and 12.7 nm (♦). The solid curves are the SDT calculated 
with the lowest two subbands included and the dashed curves 
are those calculated with only the lowest subband. 



the single subband approximation as dashed curves in 
order to reveal the contribution from the multi-subband 
effect. The figure shows that for wide QW, the SDT 
first increases with the background temperature and then 
decreases with it when the temperature rises to 150 K. 
However, for narrow QW the SDT keeps increase until 
the temperature rises to about 250 K and then decreases 
slightly. The figure also shows the contribution of the 
higher subbands: For wide QW, the population of elec- 
trons in the higher subbands is larger even in the low tem- 
perature regime. As the temperature rises, the difference 
of the SDT of multi-subband model and single-subband 
one becomes larger and larger as the contribution of the 
higher subbands becomes more important. While for 
the narrow QW, the contribution of higher subbands is 
marginal when the temperature is lower than 250 K, thus 
the SDT's from the multi- and single-subband models are 
almost the same in the regime T < 200 K. After that, 
the SDT of multi-subband model becomes smaller than 
that of single-subband one. When the temperature rises 
higher than 250 K, the SDT of multi-subband model de- 
creases with T while that of single-subband model keeps 
growing upto the room temperature. 

The different behavior of the temperature dependence 
of the SDT at different temperature regime and for differ- 
ent well width originates from the DPT. In n-type GaAs 
QW's the dominant spin dephasing mechanism is the spin 
dephasing caused by the inhomogencous broadening in- 
duced by the DPT, together with the spin conserving 
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The increase of the temperature brings 



many effects on the spin dephasing through the DPT and 
the scattering rate: Firstly, the increase of temperature 
enhances the electron-phonon and the electron-electron 
scattering, which tends to drive electrons to a more ho- 
mogeneous state, and thus reduces the spin dephasing.^" 
Nevertheless, the increase of the temperature also en- 
hances the effect from the DPT. One can see from Eq. (2) 
that the DPT in the QW are composed of the terms lin- 
ear and cubic in wave vector k = (k^, ky,0). The effects 
from the both terms increase with temperature as elec- 
trons are driven to larger wave- vector states. This tends 
to bring a faster SDT. 

For narrow QW's, the linear term is the dominate one 
even for the lowest subband as (fc^) is much larger than 
k^ in the electron density and temperature regime of our 
investigation. Although the effect of the linear term in- 
creases with temperature, the increase rate is slower than 
that of of the scattering. Therefore in the single-subband 
model the SDT of the narrow QW keeps increasing with 
temperature. However, for wide QW's, the linear term 
of the lowest subband and the cubic term of the DPT 
are comparable. When the temperature rises to 200 K, 
the contribution of cubic term becomes larger than that 
of the linear one of the lowest subband for the QW with 
a = 17.8 nm. As the effect of the cubic term increases 
much faster with temperature than the scattering, the 
SDT decreases with temperature in the regime where the 
effect of the cubic term takes over that of the linear one. 

When the contribution of the higher subbands are in- 
cluded, the increase of the temperature further enhances 
the DPT through exciting electrons to the higher sub- 
bands which brings the much larger DPT. As a result, 
for the QW with a — 17.8 nm, the SDT increases only in 
the regime T < 150 K and then decreases when the effect 
of the multi-subband is included. Similarly, the SDT of 
the 12.7 nm QW no longer increases monotonically with 
the temperature but first increases with temperature and 
then decreases near room temperature. 

We now study the effect of the applied electric field 
on the spin precession and spin dephasing. In Fig. 4, 
we plot the SDT as a function of the applied electric 
field for two temperatures and two well widths. It is 
seen from the figure that for low temperature, the SDT's 
first increase with E for both wells just as the case in 
the single-subband model.® With further increase of the 
electric field, Tg of the wide QW decreases rapidly with 
E once it is larger than 0.5 kV/cm, while that of the 
narrow QW remains growing until E is higher than 1 
kV/cm. For high temperature, the SDT decreases with 
E even in small electric field regime for the wide QW. 
While it is insensitive to the applied electric field in the 
regime E < I kV/cm for the narrow QW. 

It is understood that the electric field can affect the 
spin dephasing in two competing ways:® On one hand, 
when the electric field increases, the scattering rate in- 
creases due to the hot-electron effect. Consequently the 
electrons are driven to a more homogeneous state in k- 
space and thus the spin dephasing is reduced — Effect I; 
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FIG. 4: SDT vs. the applied electric field E at different tem- 
peratures and well widths: •, T = 120 K; ♦, T = 200 K; Sohd 
curves, a — 17.8 nm; Dashed curves, a = 12.7 nm. Inset: The 
corresponding electron temperature Te as a function of the 
electric field. 



On the other hand, the applied electric field also drives 
the electrons to higher momentum states or higher sub- 
bands by rising the electron temperature and the drifting 
velocity and thus the spin dephasing is enhanced through 
the strengthening of the DPT — Effect II. In the low tem- 
perature and low applied electric field regime, where the 
linear DPT is more important and the electron drift ve- 
locity is small, the enhancement of the DP effect is rel- 
atively smaller and the reduction of the inhomogcneous 
broadening due to the hot-electron effect is more pro- 
nounced. Therefore the SDT is raised in this regime. 
However, when the electric field becomes larger, more 
electrons are driven to high momentum states where the 
cubic DPT dominates, moreover electrons are easier to be 
excited to the higher subbands where the linear DPT is 
much larger than that of the lowest subband, the second 
effect becomes more important and the SDT decreases 
with E consequently. 

For wide QW's, electrons are much easier to be excited 
to higher subbands and the domination of the cubic DPT 
is much easier to achieve than narrow QW's. Therefore 
at low temperature the SDT of the wide QW starts to 
decrease at relatively smaller electric field (0.5 kV/cm in 
the figure). While the SDT of narrow QW keeps increas- 
ing until the electric field is 1 kV/cm, which is about two 
times as large as that of the wide well. At high temper- 
ature, as the contribution from the cubic DPT is more 
important and also more electrons are excited to high 
subbands. Effect II is more important even at low elec- 
tric field for wide QW's. As a result the SDT of wide QW 
decreases monotonically with the electric field. For the 



narrow QW in our calculation, both effects of the applied 
electric field is comparable in the regime E < 1 kV/cm, 
consequently the SDT is insensitive to the electric field 
in this regime. 

It is also interesting to note that at very high electric 
field regime {E > 2 kV/cm) , the SDT is only determined 
by the electric field and is insensitive to the well width 
and the background temperature. In this regime, elec- 
trons are excited to high energy states by the electric 
field and the electron temperature Te is much higher than 
the background temperature T as shown in the inset of 
Fig. 4. In this situation, the major scattering mechanism 
is the electron-electron Coulomb scattering and the dom- 
inant term of the DP term is the cubic one. Both rely on 
the electron temperature and are insensitive to the well 
width and the background temperature in the high elec- 
tric field regime. Nevertheless it is noted that notwith- 
standing the fact that the electric field goes to very high 
in Fig. 4, the inter-valley scattering between T-L valleys, 
which plays important role in the high-field transport, is 
not included. Once this additional scattering is incorpo- 
rated, the electron temperature will be reduced and the 
spin dephasing time may differ from what is predicted in 
Fig. 4 when E > 2 kV/cm. However a detailed study of 
this inter-valley effect is beyond the scope of this paper 
and also the spin dephasing mechanism in the L (and 
also X) valley is not full understood to date. 

In conclusion, we have studied the multi-subband effect 
on spin precession and spin dephasing in n-type GaAs 
QW's where electrons may occupy more than one sub- 
band. Our results show that due to the strong inter- 
subband scattering, electrons hop frequently among the 
subbands and thus experience an average DP effective 
magnetic field regardless which subband the electrons are 
located. Therefore the spin precession and the spin de- 
phasing of electrons in each subband are almost identi- 
cal even though the DPT's in different subbands vary 
greatly. In the system we study, the spin dephasing 
is determined by the joint effects of the inhomogcneous 
broadening induced by the DPT and the spin-conserving 
scattering such as the electron-phonon and the electron- 
electron Coulomb scattering. The scattering tends to 
randomize electrons in the momentum space and reduces 
the inhomogcneous broadening and consequently leads to 
the rise of the SDT. Whereas the enhancement of the con- 
tribution of the DPT boosts the inhomogcneous broad- 
ening and therefore reduces the SDT. When the distri- 
bution of electrons varies due to the change of the well 
width, the background temperature and the applied elec- 
tric field, the contributions from the linear and the cubic 
DPT's and the contribution from the higher subbands as 
well as the contribution from the spin conserving scatter- 
ing change accordingly and the SDT varies consequently. 

For narrow QW where electrons occupy only the low- 
est subband when the temperature is not too high, the 
dominant contribution from the DPT is the linear term. 
Therefore the SDT increases with temperature almost 
upto the room temperature as a result of the fast increase 



in the scattering rate and the relatively slower increase of 
effect from the linear DPT. While for wide QW, the con- 
tributions from the cubic DPT and the higher subbands 
are more important, the SDT only slightly increases with 
temperature in low temperature regime and then de- 
creases with it. Similar effect also appears in the applied 
electric field dependence of the SDT. In low temperature 
and low electric field regime, the SDT increases with the 
electric field as the hot-electron effect effectively reduces 
the inhomogeneous broadening of the DPT. When the 
background temperature and/or the applied electric field 
increase, the SDT decreases as the contributions from 



the higher momentum states and high subbands become 
important. Whereas for wide QW's, as the contributions 
from the cubic DPT and the high subbands are more im- 
portant, the SDT begins to decrease at smaller electric 
field than that of narrow QW's. 
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